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Abstract 

Wc present the results from 1.5D diffusion simulations of the Kippenhahn-Sehlueter prominence model 
magnetic field evolution under the influence of the ambipolar terms of Cowling resistivity. We show that 
initially the evolution is determined by the ratio of the horizontal and vertical magnetic fields, which 
gives current sheet thinning (thickening) when this ratio is large (small) and a marginal case where a new 
characteristic current sheet length scale is formed. After a timespan greater than the Cowling resistivity 
time, the current sheet thickens as a power law of t independent of the ratio of the field strengths. These 
results imply that when Cowling resistivity is included in the model, the tearing instability time scale is 
reduced by more than one order of magnitude when the ratio of the horizontal field to the vertical field is 
20% or less. These results imply that, over the course of its lifetime, the structure of the prominence can 
be significantly altered by Cowling resistivity, and in some cases will allow the tearing instability to occur. 
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1. Introduction 

Quiescent prominences are globally highly stable struc- 
tures that exist in the corona for days or even weeks. 
The characteristic density is lQ~^^gcm~^ and the char- 
acteristic temperature is ^ lO^ilT, which are two orders 
of magnitude larger and smaller, respectively, than that 
of the surrounding corona (Hirayama 1986). Using this 
value for the temperature the pressure scale height can 
be calculated to be A w SOOfcm. Quiescent prominences 
form over polarity inversion lines between areas of oppo- 
site magnetic flux (Tandberg-Hanssen 1995) and have a 
magnetic field strength between 3 ^ SOG (Leroy 1989) 
and a gas pressure of 0.06dync'm~^ (Hirayama 1986) giv- 
ing a plasma /3 « 0.1. The magnetic Reynolds number of 
a quiescent prominence is i?,„ w 10^ ^ 10^. 

Despite the global stability of the quiescent prominence, 
there are a number of observations that highlight that qui- 
escent prominences are locally highly dynamic. Engvold 
(1981) reported downflows along vertical filaments, inter- 
preted as downflows along magnetic field lines driven by 
gravity, of the order of a few kms~^ . Kubota & Uesugi 
(1986) reported flows in a dark fllament, observed on disk, 
that propagated predominately downward with a veloc- 
ity ~ Qkms^^ . Vorticies inside the prominence of ap- 
proximately lO^fcjTi X Vf'km in size with rotation rates 
■~ SOfcms"-'^ were reported by Liggett and Zirin (1984). 

Recent observations of quiescent prominences provided 
by Solar Optical Telescope (Tsuneta et al. 2008) on the 
Hinode satellite (Kosugi at al. 2007) have shown that 
on an even smaller small scale, quiescent prominences are 



highly dynamic structures undergoing constant evolution. 
Berger at al (2008) found dark upflows, of width ^SOOfcm, 
rising from the boundary between cavities that formed at 
the base of the prominence. These dark upflows would 
rise through a height of around lOA/m, before forming a 
mushroom cap profile. 

From the above, we can see that globally quiescent 
prominences are extremely stable structures, but on a 
smaller scale they are highly dynamic and cannot be con- 
sidered to be in equilibrium. One of the main problems in 
prominence physics is understanding how to capture such 
dynamics in prominence models. One possible reason for 
the formation of the localised dynamics inside the globally 
stable prominence could result from the partial ionisation 
of the prominence plasma. It has been shown that the 
ionisation fraction of quiescent prominences is around 0.3 
(Hirayama 1986). As this is not a fully ionised plasma, 
there will constantly be neutral flows inside the promi- 
nence, that would stop the formation of any equilibrium 
that was not dynamic. This has been studied observa- 
tionally by Gilbert et al. (2007), who found evidence 
of cross-fleld diffusion of neutral particles, which can be 
viewed as one mechanism to explain mass loss from a qui- 
escent prominence. Such observations imply that neutral 
dynamics are playing a role in the evolution of quiescent 
prominences. 

The impact of partially ionised plasma manifests it- 
self in a 1-fluid model in the form of Cowling resistiv- 
ity (Cowling 1957; Braginskii 1965). The impact of 
Cowling resistivity (composed of the resistive and ambipo- 
lar terms) on the magnetic field has already been studied 
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in a wide range of situations. The formation of sharp 
structures by ambipolar diffusion created by anti-parahel 
magnetic fields was studied by Brandenburg & Zweibel 
(1994). The inclusion of ambipolar diffusion was found to 
greatly enhance the growth rate for the tearing instabil- 
ity at points where the sign of the magnetic field reverses 
(Chiueh 1998). The three dimensional case is more com- 
plicated but ambipolar diffusion was shown to decrease 
the tearing timescale when the Lorentz force works to 
compress the plasma. Cowling resistivity has been dis- 
cussed as a potential driver in forming coronal nonlinear 
force free fields in the chromosphere from photospheric 
fields (Arber at al. 2009). It was shown that on a time 
scale ~ 15 mins the magnetic field would rearrange it- 
self forming a nonlinear force free field. Singh & Krishan 
(2010) explained how the ambipolar terms are effective at 
damping the propagation of Alfven-like wave modes. 

This study looks at how the prominence magnetic field, 
based upon the field used in the Kippenhahn-Schlucter 
(K-S) prominence model (Kippenhahn and Schlueter 
1957) evolves under the influence of Cowling resistiv- 
ity. The K-S prominence model is based on the simple 
idea of supporting plasma against gravity through the 
Lorentz force via magnetic tension. This model has been 
shown to be linearly stable for ideal MHD perturbations 
(Kippenhahn and Schlueter 1957; Anzer 1969). Low 
and Petrie (2005) developed a modified K-S prominence 
model, created by using the K-S model to describe one 
sheet inside the prominence and then using an array of 
sheets with varying properties to create an entire promi- 
nence. This prominence model produced constant flows 
associated with a structure in local equilibrium but glob- 
ally not in equilibrium that matched well with observa- 
tions of flows in prominences. Petrie & Low (2005) dis- 
cuss the effect of a constant diffusion on the modified K-S 
model and showed how reconnection could drive a net up- 
ward movement of magnetic field associated with a net 
downward flow of mass. 

The inclusion of Cowling resistivity and how this effects 
the current sheet profile of the K-S model are yet to be 
investigated. As the ratio of the Cowling resistivity to the 
Spitzer resistivity in quiescent prominences is rjc /rjKi 1000 
(i.e. the ambipolar terms dominate), the timescale of the 
deformation of the magnetic field would be greatly re- 
duced through the inclusion of neutral dynamics. Figure 
1 shows the K-S prominence model magnetic field con- 
figuration with a fiuid made up of neutral and charged 
particles. As the neutrals do not directly feel the mag- 
netic field, they will be able to flow out of the current 
sheet across the magnetic field (dashed arrows show di- 
rection of motion) , as the is no longer an equilibrium the 
charged particles are advected by the Lorentz force (solid 
arrows). The evolution of the magnetic field of the K-S 
model under Cowling resistivity will provide an insight 
into the nature of prominences and whether it is possible 
for them to form an equilibrium. This can then be used to 
shed some light on how to connect the large scale stabil- 
ity created in prominence models with the local dynamism 
observed in prominences. 



In §2 we present the numerical method and in §3 we 
present the results and use these to perform an estimate 
of the tearing timescale in the K-S model. In §4 we present 
a summary and a discussion outlining the physical mean- 
ing behind the results and their implications for quiescent 
prominences. 

2. Numerical Method 

The equations for the magnetic field in the K-S promi- 
nence model are as follows: 



Bz{x) :=: Bzootanh 
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where the x direction is across the prominence (a: = is the 
centre of the prominence) , the y direction (the direction of 
the current) is along the prominence and the z direction is 
the vertical direction, B^q — Bx{Q); B^oo — Bz{oo) and A 
is the pressure scale height. As this study focuses purely 
on the effect of Cowling resistivity on the magnetic field, 
so the pressure scale height does not have any physical 
meaning in this setting. To address this, we normalised 
the characteristic length scale of the current sheet to be 
equal to L, i.e. L ~ 25^.0 A/B^oo • This gives 



Bz{x) = Bzootanh ( — ] 



(3) 



The induction equation including Cowling resistivity is 
as follows (Leake and Arber 2006): 



dB 

'dt 



= V(vxB-77J|| -i]cj 



(4) 



where ?; is the Spitzer resistivity and r/c ~i] + {^n/a)'B'^, 
where ^„ and a denote the neutral fraction and the friction 
coefficient respectively. We use a 1.5D {d/dz = 0) version 
of this equation, and taking into account that there is no 
jii component due to the geometry of the problem and that 
Bx must remain constant to preserve V ■ B = 0, equation 
4 becomes 

— .-(vxB-,cj) (5) 

This work is only concerned with the effect of the diffu- 
sion, so the effect of advection of the magnetic field is 
neglected. This assumption can be justified as initially 
the system is in equilibrium, so for the initial t = tc then 
the dynamics will be dominated by this term. This is sim- 
ilar to processes observed in molecular cloud fragmenta- 
tion under Cowling resistivity (Kudoh & Basu 2008). For 
any evolution after this time, the results from this paper 
will only provide a handle on the dynamics of the system. 
Also when using typical values for quiescent prominences 
Vc/v = 1000, so it is possible to neglect 77. Therefore, this 
work only presents the solution to the following equation 



dB, 

dt 






(6) 



The Cowling resistivity profile is taken to be rjc = 
(^^/a)B^ = r;oB^ apart from in section (3.1.2) where a 
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spatially and temporally variable profile is considered. 
The time scale is normalised to the Cowling resistivity 
timcscale r^^ = L"^ /rjc, the length scale to L and the re- 
sistivity by 77oB^. 

To solve the induction equation, we use a central differ- 
ence scheme that forms part of a two step Lax-WendrofF 
scheme based on the scheme presented in Ugai (2008). 
The calculations were performed on a grid with 10^ grid 
points, unless otherwise stated, where Ax = 0.05. The 
outer boundary is a free boundary and the inner bound- 
ary a reflective symmetric boundary, in the appendix the 
case where the outer boundary is a fixed boundary is dis- 
cussed. 

3. Results 

3.1. Current Sheet Evolution 

To analyze the temporal evolution of the magnetic field, 
we perform a series of simulations on a magnetic field as 
described by equation 1 and 3, for a range of values of -B^.. 
We use a domain of size X ~[Q, 50L], over a timespan of 
- lOOrc. 

Figure 2 shows the evolution of the current sheet and 
B^{x) for (top) Bj;o/-B^oo =0.1 and (bottom) Bm/B^oo = 
0.5. The solid, dashed, dotted and three dots -I- dashed 
lines show the distribution at t = 0, 1, 10 and 50 respec- 
tively. The thick dashed line denotes the power ~ x~'^/'^ 
for the current and ~ x^^^ for the magnetic field. The 
thick horizontal line shows Bxo/ Bzoo- 

It can be clearly seen that the current sheet forms two 
distinct sections. The first where x S [Q,Lbx\, i-e. B^ > 
Bz{x), and the second where x S [L_b,y,oo], i.e. Bz{x) > 
Bx . The reason for this split in behaviour comes from the 
definition of the Cowling resistivity, 

ilc{x)^{Bl + B,{xf)r^o. (7) 

As the term B^rjQ is a constant, it will work as a uniform 
resistivity which drives current sheet thickening. Whereas 
the Bz{x)^rjQ term is a variable in x, so if the field had no 
Bx component, this would drive thinning of the current 
sheet (Brandenburg & Zwcibcl 1994). It can be shown by 
simple integration that such diffusion results in the current 
sheet forming a power law distribution oc x^"^'^ , and as a 
result the magnetic field forms a power law distribution 
where Bz{x) ex cc^/'^. As both these processes are working, 
if Bxo < Bzoo there must be some point across the current 
sheet where they are equal, forming a new characteristic 
width of the current sheet. 

For the case where Bxo/Bzoc — 0.1 (figure 2 top), cur- 
rent sheet thinning occurs and after t = tc the current 
sheet has formed the power law distribution described 
above. Throughout the evolution of the current sheet, 
this power law component can be clearly seen at all times. 
The same can be said for the magnetic field distribution. 
For the case where Bxo/Bzoo ~ 0.5 (figure 2 bottom), cur- 
rent sheet thickening and current sheet thinning behaviour 
can both be observed. The inner part of the current 
sheet approximately undergoes constant resistive thick- 
ening, maintaining a flat current profile, and the outer 



part locally appears to form a power law distribution but 
not to the same extent as the Bxa/B^oo — 0.1 case. The 
current sheet thinning behaviour docs not continue for an 
extended period of time, as shown by the later stages of 
the evolution in figure 2 bottom. Though the dynamics 
displayed by the magnetic field are not those of a constant 
resistivity process, no clear power law distribution forms. 

We found that the temporal evolution of this new char- 
acteristic width, and as a result the magnetic field, can 
be divided into two separate stages. The initial stage as 
that where approximations of the evolution based on the 
initial magnetic field distribution still hold. 
3.1.1. Initial Evolution 

The initial evolution of the magnetic field is determined 
by the ratio Bxo/Bzoo and the initial current sheet half 
width L. Inserting the definition of Cowling resistivity 
(equation 7) into the induction equation (equation 6), 
gives 



^ = -s[(«-+«=<^«™l 



(8) 



Assuming that 3x s.t. dBz{x)/dt = 0, it is possible to 
form a comparison of the two components of the diffusion 
term. 



rioB^Jy ^r]oBz{x)'^Jy 
This equation can be simplified to give 
Bx 



Bz 



' tanh 



(!) 



(9) 
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from the definition of Bz . We define the value of x found 
here {x = Lbx), as a new characteristic width of the cur- 
rent sheet, where 



Lbx = Ltanh 



B 



xO 



B, 
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B, 



(11) 



it can be seen that this lengthscale is independent of the 
resistivity profile assumed. Figure 3 shows the evolution 
of Lnew in time. It can clearly be seen that this approx- 
imation holds for a timespan of i ~ tc and for the range 
0.3 < Bxo/Bzoo < 0.7. 

The above range results from comparing the following 
terms: 

MAX {BlJy) = MAX (BzixfJy) (12) 

As Bx is constant, MAX (BlJy) = BlMAX{Jy) = Bl 
giving 



Bl = MAX (BzixfJy) ~ 0.25^2 



(13) 



which gives Bx/ Bzoo ~ 0.5. Therefore any marginal case 
will occur when this relation approximately holds. Such 
evolution is well described by the Bx/Bzoo = 0.5 case. As 
can clearly be seen in figure 2, the point at which Bz{x) = 
Bx initially, the value aX t = tq and the linear estimate 
are approximately the same value. 

As stated above, this approximation falls down in the 
case where Bl/ Bl^ < 1 and Bl/B^^ > 1. As the 
Cowling resistivity is defined as rjc = {Bl + Bz{x)'^)r]o, 
when Bl/Bl^ <^ 1 then r]c ~ Bzix^-q^. In this case we 
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find that the Hiicar stage is defined by current sheet thin- 
ning as found by Brandenburg & Zweibel (1994). This 
case is well followed by the Bx/B^oo = 0.1 case, as shown 
in figure 2. Though the prediction matches well with the 
initial value of x where B^ix) = B^, the value at < = r^ 
is completely different. The current sheet undergoes thin- 
ning, and at t = tq we have a current sheet half width 
that is 10% of the prediction. 

For the case where 1.0 < B^q/B zoo the linear approx- 
imation does not hold any more as the new characteris- 
tic current sheet half width can only be defined for the 
case where B^q < B^oo, because Va; : Bz{x) < Bzoo- The 
behaviour of the system can be understood in the same 
way as used when B^^/Bl^ ^ 1, if B^^/B^q < 1 then 
rjc ** B^riQ. This works in a similar fashion to a constant 
Spitzer resistivity, which would diffuse the current sheet 
causing thickening to occur. 
3.1.2. Later Evolution 

In the later stage, the presence of the constant horizon- 
tal field dominates the dynamics. It can be seen in figure 
3 that the behaviour of the system changes dramatically 
from the initial stage. In all cases current sheet thickening 
occurs independent of the value of B^o/ B^oo- The current 
sheet thickening is driven at a rate Lnew oz t^", which is 
analogous to current sheet thickening through uniform re- 
sistivity applied to a 6 function current sheet. It has been 
shown that the half width of a (5 function current sheet's 
temporal evolution under constant resistivity can be de- 
scribed by L = 2(77ot)^/^ (see, for example. Priest (1982)). 

The results presented so far have assumed a con- 
stant Cowling resistivity. In a quiescent prominence, the 
Cowling resistivity would vary as a function of x as the 
density and ionisation fraction vary from coronal values 
to the values at the heart of the prominence. To create 
a simple representation of this, under the assumption of 
constant gravity and hydrostatic equilibrium we have: 



B^odBzix) ^ . . 
pg = -;— — — ex Jy(X). 



(14) 



Att dx 

From this we can assume that r/ ex Jy{x) can provide an 
approximation of the prominence Cowling resistivity pro- 
file. Figure (4) shows the comparison in behaviour be- 
tween the constant resistivity and the variable resistivity 
case. As predicted the initial evolution does not depend 
greatly on the resistivity profile. The later evolution, the 
system takes a power law distribution oc i°-^. Even with 
this, it can be seen that up until t = lOrc (approximately 
the lifetime of a prominence) , the current sheet width for 
both systems is approximately the same. This allows the 
results obtained with the constant resistivity to be applied 
as general results. 

A discussion about the effect of the boundary is pre- 
sented in the appendix. 

3.2. Tearing Instability 

The formula to calculate the growth rate for the tear- 
ing instability given by FKR theory (Furth ct al. 1963) 

is (jJtear OC T^^ a~'^'^Rm* , whcrc the * symbol denotes 
that this is the effective value, i.e. calculated using the 



current sheet half width as the length scale and a = ka. 
Steinolfson & van Hoven (1983) showed numerically that 
if the constant VP approximation is not used the growth 
rate is uitear oc TAtfO^i'^Rrn* ■ Usiug the most unstable 
wavelength, a,nax oc Rm*^^'^, the fastest growth rate be- 

— 1/2 

comes uJtear OC TA*Rm* ■ If wc take Bxo/Bzoo = 0.1 for a 
current sheet half- width defined as 2BxoA./Bzoo, i-e. that 
used in the K-S model, and a global magnetic Reynolds 
number of Rm = 10^, we get a tearing time of Ttear oc IC^s. 
Using the new current sheet half width found in this 
work, Lnew = BxoLoid/Bzoo, the tearing time becomes 
Ttear = 3G0s. This rcsult, though, does not take into ac- 
count the effect of the stabilizing effect of the horizontal 
field. 

The effect of the horizontal field, i.e. a component of the 
magnetic field the is across the current sheet, was studied 
in relation to the K-S prominence model by Nishikawa 
& Sakai (1982), where it was found that for values 
Bxq/Bzoo > 0.1 the growth rate of the instability scales 
as w oc T^^ R^^\ which would give a tearing time of the or- 
der of lO^s in quiescent prominences. A similar result was 
obtained by Harrold et al. (1995) in relation to the tear- 
ing instability in the Earth's magneto-tail. Their result 
implies that as B^o/Bzoo gets larger, the transition of the 
growth rate from oc Rm* to oc R^l, happens at smaller 
values of Rm*- For example, for B^o/Bzoo = 10~^ this 
transition happens at Rm* = 10^. As such weak horizon- 
tal fields would be unrealistic in quiescent prominences, 
it is very unlikely that the tearing instability could occur 
without an outside process driving it. Sakai & Nishikawa 
(1983) showed that the tearing instability could be excited 
in such magnetic fields if driven by an outside force, in 
their case provided by short wavelengths fast mode waves. 
We believe that thinning of the current sheet caused by 
Cowling resistivity could also allow the tearing instability 
to occur. 

Using the new current sheet widths found, it will be 
possible to make an estimate for the tearing instabil- 
ity timcscale in a quiescent prominence, this estimate is 
shown in figure 5. The solid line denotes the estimate 
for the tearing timescale based on the current sheet half 
width from the K-S model, the dashed line denotes the 
tearing timescale based on the linear estimate for the cur- 
rent sheet half width under Cowling resistivity and the 
stars denote the tearing timescale for the current sheet 
half width found from the simulation results, these are 
taken at when t ~ tc. 

Figure 5 clearly shows that the inclusion of Cowling re- 
sistivity reduces the tearing timescale when B^o/Bzoc < 
0.5. The Cowling resistivity timescale is approximately 
lO^s, during which time we would expect the magnetic 
field to have restructured itself due to the flow of the neu- 
tral particles, therefore any tearing time that is longer 
than this can be viewed as unrealistic. From these results, 
when we include Cowling resistivity, we can say that for 
the case where B^/Bzoo < 0.2, it is possible that tearing 
could occur. Whereas, when wc do not include the effect of 
Cowling resistivity, it is only possible where B^/Bzoo ^ 1, 
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that is to say, unrealistically small values for the horizontal 
field component. Based on this result, we feel that the in- 
clusion of neutral dynamics could be very important when 
studying quiescent prominences. However, when applied 
to active region prominences, i.e. those with strong hor- 
izontal fields, the timescale for the instability is around 
lO^s so would not be able to occur inside the life time of 
the prominence, so such dynamics may not be important. 

4. Discussion and Summary 

One of the big problems in understanding quiescent 
prominence dynamics is explaining the highly dynamic 
features that are seen in the prominence using a model 
that has been designed to explain the global stability of 
the prominence. Until now, it has been commonplace to 
ignore the physics that are involved with partially ion- 
ized plasmas. Using values that are typical for quiescent 
prominences gives timescale for these effects of approxi- 
mately 1.5 days, which is shorter than the average lifetime 
of a quiescent prominence. Therefore the physical effects 
from partially ionized plasma could be very important for 
creating the highly dynamical features seen in quiescent 
prominences. 

In this work we have shown the impact of including the 
neutral component in a 1-fluid model, by using the addi- 
tion of the Cowling resistivity to the induction equation, 
to understand how the magnetic field in the Kippenhahn- 
Schlueter prominence model behaves. We have shown that 
the temporal evolution of the magnetic field, characterized 
by the evolution of the current sheet, which is defined 
by the value of x such that Bz{x) = B^o The evolution 
can be divided into two clear regimes, the initial evolu- 
tion which can be determined by the ratio B^o/Bzoo and 
the later evolution where the current sheet thickens with 

T nr ^1/2 

Without performing full 3D MHD simulations that in- 
clude the Cowling resistivity, it is impossible to fully de- 
scribe the dynamics that would be created inside the K-S 
model, though these relations provide a handle on how 
the model will evolve. As the plasma will flow in the 
opposite direction to the Lorentz force, for x < Lb^ the 
tension force dominates the magnetic pressure in the K-S 
model. This implies that the plasma will flow downward. 
For X > Lb^ the magnetic pressure force dominates the 
tension force, therefore there will be an outflow occurring 
along the x-direction. 

This study has allowed estimates of the tearing time 
scale to be calculated for the K-S model, giving the ratio 
Bx/Bzoo at which tearing is possible. It is clear that the 
inclusion of Cowling resistivity radically reduces the tear- 
ing time scale in the K-S model. As a result, the range of 
values of Bx/Bzoo at which we are likely to have tearing 
occur increases from ~ 0.01 to ^ 0.2. This implies that 
any numerical simulations of the tearing instability in so- 
lar prominences should include Cowling resistivity. This 
result, when applied to solar prominences, could give a 
handle on the magnetic field strength at positions where 
plasmoid ejection is observed. 



It is very important to extend this study to include fluid 
effects by solving the full MHD equations. The results in 
this study for the initial evolution of the magnetic field 
(0 < f < 1) should match well with any results from MHD 
simulations, but after that the fiuid effects should become 
important. This is because tc/ta ~ 10000. Therefore, 
once the flows have been initiated by the decoupling of the 
magnetic field from the plasma via the Cowling resistivity, 
these flows could have a signiflcant effect on the results. 
In this work, we found that the later evolution of the 
magnetic fleld for rjc = r]o was dcflned by a current sheet 
thickening with Lnew oc t^'"^ , without the inclusion of fluid 
flow this can only be viewed as a handle on the potential 
dynamics of the system. In reality, 2-fluid simulations 
with ionisation from coronal x-rays would be necessary to 
truly understand the dynamics of the system. How this 
is effected by fluid flow would be a very interesting future 
research topic. 
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Appendix. Boundary Effects 

In the main body of this work, the temporal evolution 
of the current sheet was discussed. However, for numer- 
ical simulations it is important to know the effect of the 
boundary on the system. As stated previously, a free 
boundary condition {dBz/dx = 0) was used for the cal- 
culations presented in this work. We can simply calculate 
the steady state formed with this boundary condition by 
rewriting equation 4 and neglecting the Spitzer resistivity. 
We can write 



— - = V X [(v X vc) X B] 



(Al) 



where vc = (Cn/<^)j ^ B. Therefore the effect of Cowling 
resistivity can be understood to be equivalent to advect- 
ing the magnetic field in the direction of the j x B vector, 
which in terms of the equation of motion is the Lorentz 
force. Therefore, in general, any steady state of the mag- 
netic field will be the state where j x B = 0, if the boundary 
allows such a state to be formed. 

The free boundary condition used in this study allows 
the magnetic field to relax to a j x B = state, i.e. Vx : 
Bz{x)=0 leaving a constant horizontal field. A symmetric 
boundary would allow the same process to occur. A fixed 
boundary {dBz{LBND)/dt = where Lbnd is the value 
of x at the boundary) would not allow such a state to 
be reached, because the angle of the magnetic field at 
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the boundary is fixed. Therefore, any steady state that 
is reached will be such that 3x such that j x B 7^ 0. In 
this subsection, we will investigate the steady state formed 
when using a fixed boundary. 

Figure 6 shows the temporal evolution of L^x for a 
solid boundary (Lbnd) at x = lOL & 20L. The dashed 
and solid lines correspond to the case where Bx/B^oo = 
0.4 for Lend = ^^L and 20L respectively, and the thick 
dashed and thick solid lines correspond to the case where 
Bx/Bzoo = 0.7 for Lbnd = lOL and 20L respectively. It 
is clear that once the steady state has been reached a 
characteristic current sheet half width is formed that it 
has a dependence on the value of Lbnd and B^/Bl^. 
Analysis of the results shows that the current sheet half 
width of the steady state L steady oc LbndB^J B^^. 

As shown in section 4, a steady state can be formed 
whilst the field lines remain curved. This is only ap- 
plicable for a fixed boundary, which docs not apply to 
Solar prominences as the magnetic field as the surround- 
ing corona would act like a free boundary due to its low 
density. This type of magnetic field equilibrium could be 
applicable to the emergence of flux tubes or sheets that 
rise from below the photosphere into the solar atmosphere, 
where the photosphere can act as a fixed boundary. The 
small ionisation fraction in the chromosphere means that 
if the emergence of the fiux tube is sufficiently slow, the 
structure of the magnetic field could be greatly altered to 
form an equilibrium. 
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Fig. 1. K-S magnetic field with neutral and charged particles. Sohd arrows show Lorentz force and dashed arrows show direction of 
neutral motion. 
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Fig. 2. Plot of jy (Left) and B, (Right) for B^ = 0.5 (top) and 0.1 (bottom). The solid, 

dashed, dotted and three dots + dash lines show the distribution at t = 0, 1, 10 and 50 
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Fig. 3. Plot of the evolution of the value of x = Lgx, i-c. the value of x where B^{x) = Bx, for Bx/B^oo = 0.1 (solid), 0.2 (dot 
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Fig. 5. Plot of the of tearing time scale against t(sec) for different initial values of Bxo/ B,ao- The solid line is for the original K-S 
model, the dashed line show the tearing timescale based on the linear estimate of the current sheet width under Cowling resistivity 
(see equation 11) and the stars denote the tearing timescale for the current sheet width found from simulations (see fig. 3) at t = tc 
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Fig. 6. Plot of the evolution of the value of x = Lbx, i-c. the value of x where B^{x) = Bx, for Bx/Bzoo = 0.4 
(dashed-LsjvD = lOL and solid-LsjvD = 20L) and 0.7 (thick dashed-Lgivo = lOL and thick solid-Lsivo = 20L) 



